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COMPARISON OF SEVERAL METHODS FOR OBTAINING THE TIME RESPONSE OF LINEAR 
SYSTEMS TO EITHER A UNIT IMPULSE OR ARBITRARY INPUT FROM 

FREQUENCY-RESPONSE DATA 1 

By Jambs J. Donegan and C-abl R. Hubs 


SUMMARY 

Several methods oj obtaining the time response of linear 
systems to either a unit impulse or an arbitrary input from 
frequency-response data are described and compared. 

Comparisons indicate that all the methods give good accuracy 
when applied to a second-order system; the main difference is 
the required computing time. Several of the methods when 
applied to higher order systems require excessive computing 
time in order to obtain the same degree of accuracy. The 
methods generally classified as inverse Laplace transform 
methods were found to be most effective in determining the 
response to a unit impulse from frequency-response data of 
higher order systems. 

Some discussion and examples are given of the use of the 
methods as flight-data-analysis techniques in predicting loads 
and motions of a flexible aircraft on the basis of simple calcu- 
lations when the aircraft frequency response is known. 

INTRODUCTION 

The frequency-response type of analysis used on linear 
systems has found extensive application in the field of air- 
craft stability and in the determination of overall dynamic 
characteristics of an aircraft. In fact the current trend 
to perform analysis on flight test data in frequency-response 
form appears to be gaining favor especially in the case of 
flexible aircraft. Considerable emphasis has, therefore, been 
given to methods of determining the frequency response of 
a system from transient responses. A concise r&um6 and 
comparison of methods for obtaining the frequency response 
from transient responses are presented in reference 1. 

This report briefly considers the next step in the process — 
the converting of the information contained in a frequency 
response to the time plane in the form of the response to 
a unit impulse. For certain purposes this conversion gives 
data in more useful form. The response of a linear system 
to a unjt impulse may be used in conjunction with Duhamel’s 
(superposition) integral to determine (1) the aircraft tran- 
sient response to any type of input or (2) the input required 
to cause any required aircraft transient response. 

It appears that methods of converting frequency-response 
data to transient data present a flight-data-analysis tech- 
nique which permits the prediction of aircraft motions 
and loads for a flexible aircraft without knowledge of the 

• Supersedes NAOA Technical Note 3701 by James J. Donegan and Carl R. Hues, 1968. 
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equations of motion relating the input and output. These 
methods also bypass the need for computing transfer- 
function coefficients or stability derivatives in predicting 
these loads and motions. Such predictions are important 
in anticipating the motions and loads for more hazardous 
aircraft maneuvers. These methods are also useful in 
predicting time responses of complicated linear systems 
whose frequency response is known. 

The purpose of this report is to collect and briefly compare 
a few of the methods now available for performing this 
operation. The methods are compared on the basis of 
accuracy, computing time required, and applicability of 
the method to higher order systems. Some extensions of 
these'methods are also given. With such information avail- 
able, engineers may then select the method which best 
fits their needs. 


SYMBOLS 

F{t) function of time 

H(iw) frequency response of a system 

H(s) transfer function (in terms of Laplace variable s) 

relating input and output 
time response of a system to a unit impulse 


index of summation 

transfer coefficients defined by equation (8) 
limits on summation of P-transform 
Mach number 

index of summation of P-transform 
total number of data points used in solution of 
equation (9) 

n limit of summation 

P polynomial transform operator 

Re [B'(uo)] real part of frequency responso defined by 
equation (21) 



r 


ordinate of pulse used to fit Re [H r (i<o)] 

s 


Laplace variable, c-j-uo 

t 


time, sec 

X 


output 

x{s) 

Laplace transform of a ft) 

x{t) 

response of a system to an arbitrary input 

X 

5 


amplitude ratio of frequency response 
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A 

indicates increment 

CD 0 

5 

input 

■O 

m 

time history of input to a linear system 

So 

6 

pitching velocity, radians/sec 

0 -100 

T 

dummy time variable of integration, sec 

CP 

c 

4> . • 

phase angle, deg 

ay 

4>n 

phase angle between output x and input 5 of 
frequency response; negative phase angles 

1 -200 


indicate lag 

.10 

Ci ) 

circular frequency, radians/sec 



damped natural frequency, radians/sec 


«/ 

fundamental frequency, radians/sec 

.08 

03 n 

undamped natural frequency, radians/sec 


Subscripts: 


><\co 

e 

elevator 

6 06 

i 

index of summation indicating row 

O 

3 

index of summation indicating column 

<£> 

"O 

Z> 


— — uivuvwvvw ^ uiuuiuiviiu \_/ JL *. UUVV1V11 y ^.VJL 

example, ^( 2 ) denotes polynomial transform of ^(t). 
DESCRIPTION AND DISCUSSION OF METHODS 

In this section the methods for determining the response 
to a unit impulse and the methods for obtaining the response 
to arbitrary inputs are discussed. 

The methods for obtaining the time response to a uni t, 
impulse from frequency-response data fall into several basic 
categories which may be generally classified as follows: 

(1) Inverse Laplace transform methods 

(2) Fourier method 

(3) Other methods 

These methods are briefly outlined to indicate the techniques 
involved. In order to describe the computations required 
and the accuracies obtained, each method is applied to the 
frequency-response data of figure 1 which defines a simple 
second-order system described by the transfer-function 


s 2 -f 6s+10 


( 1 ) 


It is usually intended that these methods be applied to higher 
order systems as will be shown later in the report. No 
attempt is made to repeat the development of the methods 
since this information may be obtained from the references. 

INVERSE LAPLACE TRANSFORM METHODS 

Floyd’s method. — The method developed by George F. 
Floyd and described in detail in reference 2 is referred to as 
Floyd's method. Floyd shows that the inverse Laplace 
transform h(t) of H(s) given by the integral 


h(t)== iiZ +, i I H ® et ‘ d8 

is for all positive values of time equivalent to 

h(t) =— He [Ml cos too den 

ttJ 0 


( 2 ) 


(3) 


This operation is based on the assumptions that H(s) may 
be written as the ratio of two rational polynomials in a with 
real and constant coefficients, that lira H(s)= 0, and that 


t ■<* 
!• 

.02 


0 


8 12 16 
Frequency, w, radians /sec 


20 


24 


Figure 1. — Frequency response of seoond-order system. 
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H(s) has no poles in the right half of the s-plane or on tho 
imaginary axis. The procedure for performing tho integra- 
tion required by equation (3) is to plot Re[j5T(ta>)] against w 
and then to approximate the exact shape with a series of 
straight-line segments. The straight-line approximation 
written as a sum of trapezoidal functions and equation (3) 
applied to each of the trapezoids; the resulting time functions 
due to each trapezoid are then added to obtain h(t). 

A simple illustration is shown in sketch A. The time 



function associated with a typical trapezoid (as shown by 
the cross-hatched section) may be expressed as 


, 2 , /sin ci> 3 t\ /sin A 2 t\ 


where 


lt>2= 


2 

. £C e — COj, 

A 2 -— 2 — 


.s .a 




METHODS FOB OBTAINING THE IMPULSIVE KESPONSE OF LINEAB SYSTEMS FBOM FREQUENCY-RESPONSE DATA 1025 


In general for n trapezoids the time response is given by 


w,. 2J5-* . /sin U]t\ /sin Ay£\ 

w<)= ;§ A \~i^rA-^r) ® 


In figure 2 the real part Re [H(u>)\ of the simple second- 
order system computed from figure 1 is shown plotted 
against u. The function h(t) is then easily evaluated since 

tables of are given in reference 2. In the upper part of 

37 

figure 2, two straight-line fits to the Re [S'(i-co)] are shown, 
one a five-line fit and the other an eleven-line fit. For 
clarity, however, only the points and not the connecting 
lines of the fit are shown. In the lower part of figure 2 the 
resulting responses to a unit impulse computed by Floyd’s 
method are compared with the exact response to a unit 
impulse. The accuracy of Floyd’s method depends on the 
number of lines used to fit the Re [2?(i<o)] and the location 
of the cutoff frequency. The cutoff frequency is defined as 
the maximum frequency at which the Re [Z7(-fio)] was fitted. 
It should be noted that the limits of the definite integral of 
equation (3) are 0 to °° ; however, in the practical case the 
Re [£?(■£«)] is cut off at some finite frequency. This error is 
reflected at the low values of time, especially at £= 0. For a 
given number of fines, the accuracy also depends on the 
judicious fit of the lines. 

Numerical-integration method. — In order to use automatic 
computing machines to perform the inverse Laplace trans- 
form method, the necessary operations indicated by equation 
(3) are performed in the following manner by using a numer- 
ical-integration method. For a given value of time t, the 




Figure 2. — Time response to a unit impulse obtained from Floyd's 
method compared to the exact values for the second-order system. 


product curve Re [2J(io>)] cos tu> is evaluated over the « 
range and integrated by numerical-integration techniques 
which give one point on the time history of the response to 
a unit impulse. By repeating the above computation for 
all the desired values of time, the time response can he 
obtained. 

This method was applied to the frequency-response data 
of figure 1 . The accuracy of the method is shown in figure 3. 
In the upper part of figure 3 the error due to the interval 
chosen for the computation is shown by the comparison 
between the circle and square symbols. The accuracy of 
the computation for the case of Aco=1.0, as indicated by the 
circle symbols, was not satisfactory beyond £=2.0 seconds. 
This result is to be expected since the interval Aco=1.0 was 
too large to permit numerical-integration methods to perform 
adequately the integration required by equation (3). 

Rectangular-pulse method. — A method for determining 
the time response to a unit impulse from frequency-response 
data has been given in reference 3 and has been recently 
extended in reference 4. It is referred to herein as the 
rectangular-pulse method. The method involves the use 
of tables for time-plane values equivalent to unit rectangular 



Time, /, sec 


Fiqube 3. — Time response to a unit impulse obtained from the 
numerical-integration method compared to the exact values for the 
second-order system. 
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pulses of the He [27(i«)]. The method requires fitting the 
Re [Z7(i«)] with a series of rectangular pulses or a staircase 
function so that the area under the curve is equal to the 
area of the pulse in each case as is shown in sketch B. This 



fitting although made visually should be made carefully. 
If the ordinate of pulse 1 is designated as n, of pulse 2 as r s , 
and of pulse n as r n , then the time response to a unit impulse 
corresponding to this Re [ff(ho)] is given by 

h(t)='Z^r i h J (t) (6) 

l-i 


The time functions hi(t), h„(t ) associated with each 

of the rectangular pulses shown in sketch B are tabulated in 
reference 3. These time functions are the inverse Fourier 
transforms of unit rectangular pulses of the Re [17 (ho)]. 

The method was applied to the frequency-response data of 
figure 1. The fitting of the staircase function to the 
Re [lT(ho)] is illustrated in the upper part of figure 4 and the 
accuracy of the method is demonstrated in the lower part of 
figure 4. Again the error in the computed response to the 
unit impulse at t—0 is due to cutting off the Re [S'(ho)] at 
co =20 radians/sec and not integrating to ® as indicated by 
equation (3). 

FOURIER METHOD 


Fourier response to unit impulse. — In • reference 5, the 
response to a unit step input is derived by the Fourier 
method. Differentiating this response gives an expression 
by which the response to a unit impulse may be approxi- 
mated. For a linear system characterized by its frequency 

\x\ 

response with amplitude ratio 
expression is 


and phase angle <fe}, this 


■JT fieri 




cos [(2n— l)«/f+ ( 2 ,-D U/ ] (7) 


In using equation (7) the choice of the fundamental frequency 
co/ determines the accuracy and length of the computations. 
Since the accuracy is affected by co/ it has been found by ex- 


perience that a value a>/=-jr can be used as a good first esti- 


mate, where co„ is the undamped natural frequency of the 
system. Instead of co„, the value of frequency at which the 
amplitude ratio peaks for a lightly damped system may be 
used. This value may be determined from the frequency 
response of the system. 

In table I the numerical computation of h(t) for the system 
defined by the frequency response of figure 1 is shown. In 
this computation, 14 terms are carried in the expression for 




Figure 4. — Time response to a unit impulse obtained from tho 
rectangular-pulse method compared to the exact values for the 
second-order system.. 


h(t), and more accuracy may be obtained, of course, by carry- 
ing more terms. The accuracy of the method is shown in 
figure 5 by the comparison between the computed and exact 
response to the unit impulse. 

OTHER METHODS 

Schumacher’s method. — The method of reference 6 per- 
mits the computation of transfer-function coefficients by 
assuming the shape of the transfer function relating the input 
and output and then curve fitting this relation to tho air- 
craft data in frequency-response form. Once the transfer- 
function coefficients are known, the system is completely 
specified since the response to a unit impulse input or any 
other arbitrary input may then be computed by the normal 
methods available for solving differential equations. To 
illustrate the method a sample computation is performed by 
using the system defined in figure 1. The differential equa- 
tion relating input 5 and output x for the frequency response 
shown in figure 1 is 

D x x J r KiDx J rKjX=K 3 5 ( 8 ) 

By applying the vector least-squares method of reference 0, 
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TABLE I 


where 


NUMERICAL CALCULATION ILLUSTRATING 
COMPUTATION OF THE FUNCTION FOR THE 
RESPONSE TO A UNIT IMPULSE BY THE 
FOURIER METHOD 

j^«/=|«n=0.64j 


n 

2n— 1 

(2 n — 1 )oif 


X 

8 


<t>xh deg 

i 

1 

0.64 

0. 0970 

-22 

2 

3 

1. 92 

. 0760 

-61 

3 

6 

a 20 

. 0520 

-91 

4 

7 

+ 48 

. 0350 

-111 

6 

9 

6. 76 

. 0240 

-124 

6 

11 

7. 04 

. 0170 

-133 

7 

13 

a 32 

.0130 

-140 

8 

16 

9. 60 

.0100 

-145 

9 

17 

10. 88 

. 0079 

-149 

10 

19 

12. 16 

.0064 

-162 

11 

21 

13.44 

. 0053 

-155 

12 

23 

1+72 

.0044 

-157 

13 

26 

16. 00 

. 0038 

-159 

14 

27 

17. 28 

. 0033 

-160 


[0.097 cos (0.64i- 22) +0.076 cos (1.924—61) -f- 

0.062 cos (3.24-91) +0.036 cos (+484-111) + 

0.024 oos (6.764-124) +0.017 cos (7.044-133) + 

0.013 cos (8.324- 140) +0.010 cos (9.64-146) + 

0.0079 oos (10.884-149) +0.0064 cos (12.164-162) + 
0.0063 cos (13.444-166) +0.0044 cos (1+724-157) + 
0.0038 cos (164-169) +0.0033 cos (17.284-160)] 


the following set of simultaneous equations is obtained: 



These equations are then solved simultaneously for the 
transfer coefficients K 3 , K u and K 2 . From the form of 
equation (8), it is seen that the response to a unit impulse 
input is given by the equation 


jr 

hit) = — - e -31 sin a d t 


( 10 ) 


and 




A sample computation is shown in table II. Table IE 
demonstrates the computational steps involved in solving 
for the transfer coefficients K u K 2 , and K 3 for the system 



Figtjbd 6. — Time response to a unit impulse obtained from the 
Fourier method compared to the exact values for the second-order 
system. 


TABLE n 

NUMERICAL CALCULATIONS ILLUSTRATING 
SCHUMACHER’S METHOD 


radians/sec 

\x\ 

U| 

deg 

A X 

- cos (pxs 

t> X ' J 

£,= - sin <£xj 

0 

a iooo 

0 

0. 10000 

* 0 

1 

.0924 

-33.7 

. 07687 

-. 05127 

2 

. 0746 

— 63.'4 

. 03336 

-. 06661 

3 

. 0555 

-86. 8 

. 00310 

-. 05541 

4 

. 0404 

-10+ 0 

-. 00977 

-. 03920 

5 

. 0298 

-116. 6 

-. 01334 

-. 02665 

6 

. 0225 

-125. 8 

-. 01316 

-. 01825 

7 

. 0174 

-132. 9 

-. 01184 

-. 01275 

8 

. 0139 

-13R 4 

-. 01039 

-. 00923 


9. 0 

-0. 91336 

-a 15483 


'a,' 


i. 97049' 

—0. 91336 

a 15220 

0 



► 

0 

-0. 16483 

0 

0. 03069 


K, 


a 15220 

Vi ✓ 


Aj=l. 0002 


Ai =6. 0023 
£j=10. 0053 
bj= 0. 992 
a=3. 0011 

h ( 0 = 1. 0010e-»-°°u' sin 0. 9994 
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defined in figure 1. These transfer coefficients are then sub- 
stituted into equation (10), and a plot of this function is 
compared with the exact response of the system to a unit 
impulse input in figure 6. 

This method gave additional information when applied to 
the second-order system since the transfer-function coeffi- 
cients were deter min ed in the process of the computation. 
A reasonable amount of computing time was required, and 
good accuracy was obtained. The method, however, 
requires previous knowledge of the form of the tra n sform 
function relating the input and output. Use of the curve- 
fitting method on the higher order system did not prove too 
effective in determining all the parameters of the flexible 
system. If the form of the transfer function of the short 
period is assumed and fitted to that portion of the aircraft 
frequency response associated with the short period then the 
method is very effective in determining the short-period 
transfer-function coefficients. 

The P-transform method. — The P-transform method, as 
described in appendix A, differs from the other methods 
presented in this report in that the time response to a unit 
impulse can be determined directly from a known response 
to a known arbitrary input. The computation is carried out 
entirely in the time domain, bypassing the frequency plane 
entirely, and does not require knowledge of the transfer 
function relating the input and output. The method also 
represents a simple procedure for using the response to a 
unit impulse and determining the response to a given 
arbitrary input. 

An example of the method is shown in figure 7 in which 
the assumed input 5(f) and output x{t) are given from which 
the response to the unit impulse h(t) is computed. A 
comparison of the computed response to a unit impulse and 
the exact response to a unit impulse is shown in the lower 
part of figure 7. 

The z-transform method. — The z-transform method as 
described in appendix A, when applied to the continuous 
linear systems, is very similar to the P-transform method. 
In fact, for a given time interval, the z-transform gives the 
same results as the P-transform for these systems. The 
computation obtained by using the z-transform is so similar 



Figure 6. — Time response to a unit impulse obtained from Schu- 
macher’s method compared to the exact values for the second-order 
system. 



Figure 7. — Time response to a unit impulse obtained from the 
P-transform method* compared to the exact values for the second- 
order system. 


to that obtained by using the P-transform that it is not 
repeated herein. 

RESPONSE TO ARBITRARY INPUTS 

The procedure for obtaining the response to arbitrary 
inputs when the response to the unit impulse is known is 
sim ply a matter of applying Duhamel's integral to the 
response, to the unit impulse and time history of the arbitrary 
input. A numerical method of performing this operation is 
shown in appendix A of reference 3. 

The response to a unit step input may be obtained by 
numerically integrating the response to a unit impulse by 
using the integrating matrix given in reference 3. In like 
manner, the response to a ramp input may be obtained by 
numerically integrating the response to a step function. 
The response to a triangular input is obtained by super- 
position of the responses to various ramp inputs translated 
along the time scale. 

The response to an arbitrary input may be obtained from 
the response to a unit impulse by the P-transform method 
by use of equation (A4) of the appendix. 
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For directly dete rminin g the time response to an aribtrary 
input from frequency-response data, the Fourier method is 
perhaps the best known. The general method is indicated 
in references 7 and 8. In general, an input which can be 
expressed as a Fourier series can be represented as 


a(f)=-Ao+ib CUXsin (11) 

B=1 ’ ’ 


The response to the input given by equation (11) of a linear 
system whoso amplitude ratio is - and phase angle is <f>xs 
can bo written as , 


x(t) — O na 



Sin [TiMft + &, B/ + (4>xs)n] 


( 12 ) 


For a few specific shapes of inputs the values of B a , C nUf! and 
of equation (12) are known. Two wave forms fre- 
quently used are the square wave and triangular wave and 
examples of these are given. 

Response to a square wave. — As indicated in reference 5 
the Fourier series for a unit square-wave input S(t) may be 
written as 

5(0=4+^ g 2 ~l 8in (2n-D«** (13) 


and the response to this unit step input of a linear system 


of which the amplitude ratio is 
4>xt is given by 


and the phaso angle is 




+ ? f; 

=o ir n =i 


ain [(2n— l)«/t+(^j)(*»-u «/] 


(14) 

The accuracy of the computation depends on the choice of 
the optimum fundamental square-wave frequency o)/ and, as 

co 

before, a suitable value is usually &>/=-w where a>, is the 

o 

lowest undamped natural frequency of the system. 

Response to a triangular- wave input. — The specific form 
of the input used in determining the time response to a 
triangular-wave input is shown in sketch C where T f is the 



total period of the input, T { is the length of the base in 
seconds, and 5 is the maximum value of the input. 

Without giving the analytical development, representation 
of the Fourier series for this triangular-wave input is given by 


where 


5(0=SC , 2*-i sin [(2» — l)«r*H-0a»-i] 


radians/sec 


(15) 

(16) 


Csa-l- 


45 


and 


**(2»-l)*( ^ 


c<)s (2n ^GD] ( - i7) 


^ b _!=180 [J- (2n — 1) (J)] (18) 


The response to this triangular-wave input of a linear system 

1*1 


of which the amplitude is 


and the phase angle is ^ is then 


*(0=Sft»-i 

«*=•! 


(2*-l) <*, 


sin [(2ft— 1 ) 0 ^ -f- 


0<2n- 1) «/+" (<£xj) (2n- 1) (19) 


Here again a suitable first approximation is Uf= 



For a higher order system such as a linear system with 
several structural modes (as indicated by peaks in the 
amplitude-ratio curve), the choice of a fundamental fre- 
quency Uf becomes difficult inasmuch as odd multiples of 
cc f must give the natural frequency &>„ and the frequency of 
each of the higher structural modes. 

If an input is nonperiodic and cannot be expanded into a 
Fourier series directly, the Fourier transform of the input is 
used. The tables of reference 3 may be used to determine 
numerically the Fourier transform of an arbitrary input. 
The procedure for computing the response of a linear system 
to such an arbitrary input is briefly described as follows. 
The Fourier transform of the output is formed by multiplying 
the Fourier transform of the arbitrary input by the fre- 
quency response of the system, which is the ratio of the 
Fourier transforms of the output and input. The inverse 
Fourier transform of this product is then taken which is the 
time response of the system to the arbitrary input. The 
tables of reference 3 may also be used to perform the opera- 
tion of the inverse Fourier transform. 


COMPARISON AND USE OF METHODS 
COMPARISON OF METHODS 

A comparison of the methods as applied to the second- 
order system, defined by the frequency response shown in 
figure 1, may be made by noting the differences between the 
exact response to a unit impulse and the response computed 
by each of the methods shown in figures 2 to 7. From these 
-figures there appears to be little difference between the 
accuracy of any one method over the others, and in each 
case greater accuracy may be obtained at the cost of more 
computation. For the comparisons shown, however, the 
computing time required by the rectangular-pulse method 
was significantly less than the time required by the other 
methods. The Fourier response to a unit impulse required 
the most time. A listing of the methods in the order of the 
computing time required for the comparisons shown in 
figures 2 to 7 is given as follows: 

Rectangular pulse (refs. 3 and 4) 

Schumacher’s (ref. 6) 

Floyd’s (ref. 2) 

Numerical integration 
Fourier reponse to a unit impulse 
(All computations were performed on a desk-type computer.) 
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A more severe test of the methods occurs when they are 
applied to linear systems of higher order. In order to 
demonstrate an application of this type, a further com- 
parison of the methods was made by applying them to the 
system defined by the transfer function 


x, .1 100 225 

S W s 2 +6s+10 s 2 -}- 0.4s-}- 100 s i +0.2s+225 


( 20 ) 


The frequency response of this system is shown in figure 8. 
The real part of the frequency response shown in figure 9 
was computed from figure 8- by using the relation 


K.e[£f (•&!>)]= 


x 

S 


cos <t> xi 


( 21 ) 


In figure 10 a comparison of the response to a unit impulse 
computed by four of the methods is shown. For Floyd’s 
method a 27-line fit to the B.e[ff(uo)] shown in figure 9 
was used. For the rectangular-pulse method, the staircase 
fit to the Re[I7(ioj)] was made with the interval Aco=l 


.15 

.10 

.05 


^O 



-.10 


-.15 
-.20 
-.25 

Frequency, w, radians /sec 

Figure 9. — Real part of the frequency response of the highor ordor 

system. 




Figure 8. — Frequency response of higher order system. 


x , . 1 100 225 

S W s’+Os+lO s»+0.4s+100 s>+0.2s+225‘ 



Figure 10. — Time response to a unit impulse obtained by four methods 
compared to the exact values for the higher order system. 


radian/sec, and for the numerical-integration method, an 
interval Aw=0.5 radian/sec was used. The response of 
the linear system defined by equation (20) to a triangular 
input shown in figure 7 of reference 3 was used to compute 
the points shown in figure 10 for the P-transform method. 
A time interval of At=0.1 second was used. The Fourier 
response to a unit-impulse method was found to be im- 
practical in this case because of the large amount of com- 
putation required to obtain any accuracy. The accuracies 
of each of the four methods illustrated appear to be 
equivalent. 

USE OF THE RESPONSE TO A UNIT IMPULSE AS A FUGHT-DATA-ANALYSIS 

TECHNIQUE 

As a further comparison of the methods applied to a higher 
order system, a typical longitudinal maneuver (M= 0.82) 
for a flexible swept-wing airplane has been analyzed. In 
this case the output was the pitching-velocity response 
6 at the center of gravity of the airplane and the input was 
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the elevator angle 8,. The analysis of this maneuver serves 
to demonstrate the use of the method of obtaining the time 
response of linear systems to a unit impulse from frequency- 
response data as a flight-data-analysis technique. The 
Fourier integrals of the output 6 and input S„ were evaluated 
by using automatic computing equipment and the methods 
of integration of product curves were used. The frequency 
response obtained by dividing the Fourier integral of the 
output by the Fourier integral of the input is shown in 


figure 11. The Re 


"1 

— (US) computed from the frequency 

-°e J 


response of figure 11 is shown in figure 12. In figure 13 a 
comparison of the response to a unit impulse computed by 
three of the methods is shown. For Floyd’s method a 


37-lino fit to the Re 



shown in figure 12 was used and 


hi(t) was computed at enough values of time to define 
adequately its shape. For the rectangular-pulse method, 


the staircase fit to the Re 




was made with the interval 


Aw=l radian/sec and, for the numerical-integration method, 
an interval of A<o=0.5 radian/sec was used. In the case 
of the numerical method, points are shown in figure 13 only 
at time intervals of 0.1 second in order to compare accuracies. 
Several attempts wore 'made to compute the response to 



Figure 11. — A typical frequency response obtained from^flight data 
of a flexible airplane at M=0.82. 



Figurb 12. — Real part of the typical frequency response obtained from 
flight data of a flexible airplane at Af<=0.82. 



Figure 13. — Time response to a unit impulse obtained by three 
methods for the typical frequency response obtained from flight 
data of a flexible airplane at Af=0.82. 


the unit impulse by the P-transform method using the 
original elevator input 8,(i) and output 6(f) shown in figure 
14. The computations were made at time intervals ra n ging 
from 0.02 second ^Af^O.l second but did not yield satis- 
factory accuracy. It appears that this inaccuracy is 
primarily due to the sensitivity of the method to small errors 
in the first few terms of the response or input. 

Time histories of the original elevator input 5, and pitching- 
velocity response 6 are shown in figure 14. The elevator 
input 8, has also been used as a forcing function with the 
response to a unit impulse shown in figure 13 to compute a 
time history of 6. The Duhamel method outlined in 
appendix A of reference 3 with a time interval of 0.05 second 
was used for these calculations. A comparison of this 
computed 6 response with the original 6 response from a 
flight test is shown, in the lower part of figure 14. This 
comparison shows the amount of error involved in the total 
computation procedure (transferring the data in time- 
history form to frequency-response form, then to the response 
to a unit impulse, and finall y to the response to an arbitrary 
input by means of DuhameFs integral). 
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Figure 14. — Comparison with the actual flight record of the calculated 
pitching-velocity response using the same inputoutput data as was 
used to obtain the typical frequency response of the flexible airplane 
at 11/= 0.82. 

Another flight-test maneuver at M=0.80 was selected 
with the other conditions approximately the same as the 
previous data in order to see how well the response to a unit 
impulse computed from one maneuver could be used to 
predict the time response to an arbitrary input from a 
different maneuver. The time histories of this maneuver 
are shown in figure 15.' Also shown in the lower part of 
figure 15 is the 6(f) response calculated by application of 
Duhamel’s integral to the response to a unit impulse given 
in figure 13, and the elevator motion for the maneuver is 
shown in the upper part of figure 15. A comparison of the 
computed and measured 6 response in the lower part of 
figure 15 gives some indication of -how well the motions of 
an aircraft can be predicted by a detailed analysis of a 
single maneuver. 

CONCLUDING REMARKS 

Several methods were compared for obtaining the time 
response of linear systems to either a unit impulse or arbitrary 
input from frequency-response data. The methods were 
compared on the basis of accuracy, computing time required, 
and applicability to higher order linear systems. 

The application of each of the methods to a simple second- 
order system indicated little difference between the accuracy 
of one method over the others, and, in general, it would be 
expected that greater accuracy could be obtained for each 



Figure 15. — Comparison of the predicted pitohing-velooity response 
with the actual flight record obtained from the flexible airplane at 
JVf=0.80. 

of the methods at the cost of more computing time. 

For higher order systems the three methods generally 
classified as inverse Laplace transform methods were most 
effective. They gave good results for a moderate amount 
of computation. All of these methods are based on the 
evaluation of the form of the inverse Laplace transform 
equation for positive values of time. 

The methods generally classified as Fourier methods gave 
good accuracy when applied to the second-order system. 
The accuracy of these methods was found to depend on the 
choice of a fundamental frequency «y. For simple systems 

a suitable value was found to be w /~y where is the lowest 

natural frequency. For higher order systems the choice of 
becomes more difficult. The Fourier response to a unit- 
impulse method was found to be impractical in the case 
of higher order systems because of the large amount of 
computation required to obtain accuracy. 

Schumacher’s method gave additional information when 
applied to the second-order systom since the transfer- 
function coefficients were determined in the process of the 
computation. This method required a reasonable amount 
of computing time and gave good accuracy. The method, 
however, requires previous knowledge of the form of tho 
transfer function relating the input and output. Use of the 
method on higher order systems did not prove too effective 
in determining all the parameters of the flexible system. A 
valuable use of the method was found in fitting the known 
short-period transfer function to only the short-period por- 
tion of the flexible-system frequency response to determine 
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the short-period transfer coefficients. 

The P-transform method is different from the other 
methods presented in that the time response to a unit 
impulse can be determined directly from a known response 
to a known arbitrary input. The computation is carried out 
entirely in the time domain and bypasses the frequency 
plane entirely. The method also represents a simple proce- 
dure for using the response to a unit impulse and determining 
the response to a given arbitrary input. The method, 
howover, when applied to the flight data of this report did 
not yield satisfactory accuracy. It appears this is primarily 
due to the sensitivity of the method to errors in the first few 
terms of the response or input. 

For all the examples computed in this report the rectan- 


gular-pulse method of NACA Technical Note 3598 required 
less computing time, for the same degree of accuracy, than 
any of the other methods used to obtain the time response 
to a unit impulse from frequency-response data. 

When an aircraft frequency response is known, it appears 
that these methods may be used as a flight-data-analysis 
technique which permits prediction of aircraft motions and 
loads without knowledge of the equations of motion relating 
the inputs and outputs for a flexible aircraft. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 16, 1956. 


APPENDIX A 


P-TRANSFORM AND z-TRANSFORM METHODS 


THE F-TRANSFORM METHOD 

Recently a new general operational calculus has been 
introduced, the elements of which were given by Tustin in 
reference 9, and the mathematics was formalized by Bubb in 
reference 10. In this calculus, which appears well adapted to 
the numerical analysis and synthesis of linear systems, a 
polynomial transform or P-transform of a function F(t) is 
defined by the equation 

P[F(t)]=P(x) 

— S M)x” 

m =—K 

= At)zr (Al) 

where F(m At) is the ordinate of F(t) at an integral multiple 
of a time interval At. It may be noted that equation (Al) 
may be obtained from the Laplace transform by setting 
t=m At and x=e~ ,A ‘, replacing the integral sign by a summa- 
tion sign, and summing over the limits ra= 0 to m= co . The 
inverse P-transform is given by 

F(t)=P~ 1 [P(x)] (A2) 

The superposition (convolution or Faltung) integral, which 
is also known as Duhamel’s integral, relates the input 5(f), 
output x(i), and response to a unit impulse h(t) (also called 
a memory function) of a linear system and is given by 

x(t) = f S(r)h(t — r)(ir (A3) 

t/ — co 

Bubb (ref. 10) shows that the P-transform of equation (A3) is 
x(x)=AtS(x)R(x) (A4) 

In the synthesis problem the input S(t) and the output x(t) 
are given and the response to the unit impulse h(t) is to be 
calculated. This calculation is performed by forming the 
P-transforms of x(x) and S(x) and dividing by ordinary 
polynomial division, x(x) by S(x), to get the P-transform of 


the response to the unit impulse 


fv(x) 


_ 1 x(x) 

At S(x ) 

= S 

m^—K 


o 


(A5) 


For practical cases, however, the summation is made over a 
finite range of t values. The inverse P-transform to h{x) 
then gives h(t), the time response to the unit impulse. An 
example of this operation is shown in table IH; the known 
input S(t) and output x(t) are shown in the upper part of 
figure 7, and the computed response to the unit impulse is 
compared with the analytic solution in the lower part of 
figure 7. 

This operational calculus is also well adapted to the solu- 
tion of the analysis problem in which the system response 
to a unit impulse h(t) and the input function S(t ) are known 
and the calculation of the system output x(t) is desired. The 
P-transforms of h(t) and 5(f) are formed and multiplied 
together by ordinary polynomial multiplication as indicated 
in equation (A4) . Since this operation is just the inverse of 
the operation shown in table HI, an illustrative examplo is 
not shown. This method has been found to be a simple and 
rapid means of applying the Duhamel process. 

The value of this operational calculus lies in the fact that 
all computations remain in the time domain and no transla- 
tion to the frequency plane is required. Also, only simple 
direct arithmetical procedures are required for solving 
practical problems. 

THE r-TRANSFORM METHOD 

The z-transformation developed originally for the analysis 
and synthesis of sampled-data systems is also applicable to 
numerical solutions of continuous linear systems, as shown 
in references 11 and 12. The z-transform is defined by the 
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TABLE HI 

NUMERICAL CALCULATION ILLUSTRATING THE P-TRANSFORM METHOD 


t , sec 

Time histories 

P-transforms 

Inverse 

P-transform 

Input,- !(i) 

Output, x(l) 

Input, a(s) 

Output, s(s) 

Response to a unit 
impulse, h(x) 

Response to a unit 
impulse, h ( 0 

0 

0 

0 

0 

0 

a 014 

0. 014 

. 1 

. 1 

. 00014 

. la: 

. 000141 

. 072s 

. 072 

. 2 

. 2 

. 00100 

. 2 x» 

. 00100s 8 

. 106s 7 

. 100 

. 3 

. 3 

. 00292 

. as 5 

. 00292s 8 

. 119s* 

. 119 

.4 

' .4 

. 00603 

. 4s 4 

. 00603s 4 

. 117s 4 

. 117 

.6 

.6 

. 01031 

. 51 s 

. 01031s 5 

. 106s 5 

. 106 

.6 

.4 

. 01537 

. 4x° 

. 01537s 5 

. 095s 6 

. 095 

.7 

.3 

. 01994 

• 3s 7 

. 01994s 7 

. 077s 7 

. 077 

.8 

. 2 

. 02316 

. 2s« 

. 02316s 8 

. 066s 8 

. 066 

.9 

. 1 

. 02466 

. 1 xf> 

. 02466z» 

. 053s 8 

. 053 

1. 0 

0 

. 02435 

0 

. 02435s 10 

. 044s 10 

. 044 

1. 1 

0 

. 02250 

0 

. 0225 Qs u 

. 030s 11 

. 030 

1. 2 

0 

. 01977 

0 

. 01977s“ 

. 026s“ 

. 026 

1.3 

0 

. 01681 

0 

. 016813“ 

. 02ls“ 

. 021 

1.4 

0 

. 01393 

0 

. 01393s 14 

. 015s 14 

. 016 

1. 5 

0 

. 01131 

0 

. 01131s 15 

. 016s 15 

. 016 

1. 6 

0 

. 00903 

0 

. 00903s 15 

. 006s 15 

. 006 

1. 7 

0 

. 00716 

0 

. 00716s 17 

-. 004s 17 

-. 004 

1. 8 

0 

. 00552 

0 

. 00552s 18 

. 012s 18 

. 012 


S(x) 

_ 1 /0+0.00014s+0.00100a*+0.002923 8 + 0.006033 4 + 0.01031s 5 + . . . + 0.00652s 18 . 
0.1 \ 0+0.1a:+0.2x’+0.3i , +0.4^+ . . . + 0.13 8 

A(s) =0.014+ 0.072s+0.106s s + 0.119s 5 -!- 0.117s 4 + . . . + 0.012s 18 


equation 

z[F(t)]=F*(z) 

=1hf(nAt)z -• (A6) 

n— 0 

Comparison of equation (A6) and (Al) indicates the similar- 
ity of the two transforms. The 2 -transformation may be 
obtained from the Laplace transformation by setting 
t=nM and 2 =e* a, ) replacing the integral sign by a summa- 
tion sign, and summing over the limits 7i=0 to n= <». The 


numerical computation applying the 2 -transform to tho 
continuous systems used in this report is so similar to that 
in which the P-transform is used that it is not repeated. 

It should be noted that in reference 13 Salzor defines tho 
symbol z by 

2=e"* A ‘ 

This definition makes the 2 -transform, when applied to 
continuous linear systems the same as tho P-transform. 
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